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Beschreibung 

Verfahren zur Bilderstellung fur einen mehrzeiligen Spiral- 
Computertomographen mit 3D-'Rekonstruktion und Computer- Tomo- 
graphie-Gerat zur Durchfuhrung dieses Verfahrens 

Die Erfindung betrifft ein Verfahren zur Bilderstellung fur 
einen mehrzeiligen Spiral-Computertomographen mit 3D-Riickpro- 
jektion, wobei ein Untersuchungsobj ekt durch ein rotierendes 
und sich in Richtung der Rotationsachse bewegendes Strahlen- 
biindel bezuglich seines Absorptionsverhaltens spiralformig 
abgetastet und die gemessenen Absorptionsdaten durch einen 
flachig ausgebildeten Detektor gesammelt werden. Die gemesse- 
nen Absorptionsdaten werden auf einen virtuellen Detektor, 
vorzugsweise im Drehzentrum der Rotation, projiziert und ge- 
filtert, und zur Rekonstruktion des Absorpt ionswertes mindes- 
tens eines Voxels werden die gemessenen und gefilterten Daten 
verwendet, die durch Strahlen erzeugt worden sind, welche 
dieses mindestens eine Voxel durchdringen . Aufterdem betrifft 
die Erfindung ein CT-Gerat zur Durchfuhrung dieses Verfahrens 
mit entsprechenden Mitteln zur Bilderstellung. 
Des weiteren betrifft die Erfindung ein CT-Gerat zur Durch- 
fuhrung dieses Verfahrens mit entsprechenden Mitteln zur Bil- 
derstellung . 

Es ist bekannt, in der Computertomographie Mehrzeilendetekto- 
ren und Strahlenbundel mit Kegelstrahl-Geometrie zur Abtas- 
tung von Untersuchungsobj ekt en, insbesondere von Patienten, 
zu verwenden und eine 3D-Bildrekonstruktion durchzuf uhren . 
Hierbei werden integrale Absorptionsdaten von Strahlen, die 
das Untersuchungsobj ekt durchdringen mit flachig ausgebilde- 
ten mehrzeiligen Detektoren gewonnen und gesammelt, wobei die 
Detektoren mitsamt einer Strahlungsquelle spiralformig um ein 
Untersuchungsobj ekt bewegt werden. Anschlieftend werden die 
Messdaten einer geeigneten Filterung unterzogen und ausge- 
wahlte gefilterte Daten ruckpro j iziert um Schichtauf nahmen 
des Untersuchungsobj ektes zu erzeugen. Unter einer 3D-Bildre- 



200210626 ^ 3 " 



2 

konstruktion wird dabei verstanden, dass zur RiAckpro j ektion 
jeweils einzelne dreidimensionale Volumenelemente (=Voxel) 
betrachtet werden und jeweils die Daten verwendet werden, die 
durch Strahlen erzeugt wurden die dieses Voxel durchdrungen 
haben . 

Bei dieser Vorgehensweise ergibt sich das Problem, dass aus 
den gesammelten Daten, die zur einfacheren Berechnung zu- 
nachst auf einen virtuellen Detektor - der auf einer Ebene 
liegt, welche die Drehachse von realem Detektor und Strahlen- 
quelle schneidet - projiziert werden, zweidimensional vorlie- 
gen und die Filterung wegen des erzeugten Rechenauf wandes nur 
in einer Dimension stattfinden soli, keine Vorzugsrichtung 
ersichtlich ist, in welche die Filterung stattfinden soli. 
Wird eine ungiinstige Filterrichtung gewahlt, so geht dies mit 
Einbufien bei der Bildqualitat her. 

Die Aufgabe der Erfindung liegt also darin, ein 3D-Rekon- 
struktionsverf ahren zu finden, bei dem die Richtung der Fil- 
terung optimal zur Erreichung einer hohen Bildqualitat ausge- 
wahlt ist. 

Diese Aufgabe wird durch die Merkmale der unabhangigen Pa- 
tentanspruche gelost. Vorteilhafte Weiterbildungen der Erfin- 
dung sind Gegenstand untergeordneter Anspruche. 

Betrachtet man das Bild eines virtuellen Detektors mit darauf 
abgebildeten Projektionen von Segmentebenen, wie sie nachste- 
hend definiert werden und wie es beispielsweise in der deut- 
schen Pat entanmel dung mit Aktenzeichen 10106398.9 (internes 
Az. 00E8563 / 2001P02305) beschrieben ist, deren Offenba- 
rungsgehalt vollinhaltlich in diese Schrift ubernommen wird, 
so zeigt sich, dass sich uber den virtuellen Detektor hinweg 
die Richtung der Projektionen verandert. Erf indungsgemaii wird 
die Filterung der Daten in Richtung dieser Projektionen auf 
dem virtuellen Detektor durchgef uhrt . Betrachtet man weiter- 
hin den Weg der Proj ektion eines Voxels auf dem virtuellen 
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Detektor bei der Rotation, so verlauft dieses projizierte Vo- 
xel durch Gebiete unterschiedlicher Ausrichtung der Segment- 
ebenenpro j ektionen . Dabei ist es besonders gunstig die Filte- 
rung der Daten jeweils entlang der Ausrichtung der Segment - 
ebenenpro j ektionen vorzunehmen, so dass pro j ektionswinkelab- 
hangig unterschiedliche Filterrichtungen entstehen . 

Betrachtet man eine 3D-Rekonstruktion, wie sie in der o.g. 

deutschen Patentanmeldung mit Aktenzeichen 10106398.9 be- 

schrieben ist, so werden in einem Spiralsegment der Lange L s 

< n+2*p max ( = 180° zuziiglich zweimal den Facherwinkel p max des 

Detektors) , beziehungsweise < 180° nach vorhergehender Paral- 

lelsort ierung, N t nt mehrfach geneigte Segmentebenen in den 

Spiralgang eingepasst . Die maximale Neigung der Segmentebenen 

R k (l<k< N t ii t ) bestimmt sich dabei aus der Bedingung, dass an 

den beiden Enden des Spiralsegments innerhalb des Messfeldes 

Strahlen fur die Ebene vorhanden sein mussen. Beispielsweise 

gilt dann fur L s =240° (mit p max = 30°) und fur L s =180° nach Pa- 
rallelsortierung N alt = |~2 • M I p\ mit M der Anzahl der Detektor- 

zeilen und p der Steigung (=pitch) der Spirale. Erf indungsge- 
mali definieren die im Spiralsegment im Abstand Aoc = L s /M be- 
stimmten Schnittlinien der Segmentebenen im virtuellen Detek- 
tor im Drehzentrum die Richtung der Filterung der Daten. 

Wird fur k=l beziehungsweise k= N t iit das Spiralsegment links- 
seitig und rechtsseitig erweitert, so wird der virtuelle De- 
tektor vollstandig abgedeckt- Auf diese Weise gelangt man fur 
jede Segmentebene Rk (l<k< N t iit) zu einem separat gefilterten 
Detektorbild D k . Diese Filteroperation geschieht in jedem der 
proj ektionsweise gemessenen Detektorbilder , wobei die Detek- 
torbilder Dk in der Riickpro j ektion in unterschiedlichen Teil- 
segmenten des Rekonstruktionsintervalls zustandig sind. 

Zur 3D-Riickpro j ektion wird ein Spiralsegment Ii der Lange 1^ = 
[-(Xmax/ oc max ] mit otmax =M • 7i/p in N t ii t uberlappende Teilsegmente 
l! k (l<k< Ntiit) der Lange L s eingeteilt, wobei deren Schwer- 
punkte sich maximal um L s unterscheiden und aquidistant ein- 
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geteilt werden. Das Segment Ii wird mittig in das Rekon- 
struktionssegment I R eingepasst . Hieraus ergibt sich fiir die 
Subsegmente I R k (l<k< N t ii t ) : 
I R k =I / k ;l<k<N tilt 

5 I R l = // <U [~ CC V max ,-(X max} 

r mit _ j mit { v \ 

1 R — l j v-^ )GC max , OC max J / 

so dass im Rekonstruktionssegment I R k mit 1 < k<N t iit das zu ei- 
nem Bildvoxel gehorige Pro j ektionsdatum im Detektorbild D k 
proj ektionsweise ermittelt wird. 

10 

Zur Veranschaulichung dieses Erf indungsgedankens ist es vor- 
jK^lfc teilhaft, das Bild eines Voxels im virtuellen Detektor in ei- 
nem Rekonstruktionssegment I R zu verfolgen. Betrachtet man 
zunachst das Spiralsegment Ir 1 und den Proj ektionswinkel a = 

15 -Omax, so sind in diesem Spiralsegment fur die Ruckpro j ektion 
proj ektionsweise die Detektorbilder D ± zustandig. Alle Voxel, 
die im Pro j ektionswinkel a = - a^x in der erf indungsgemaft 
gewahlten Filterrichtung im virtuellen Detektor abgebildet 
werden, bilden eine gekriimmte Flache Hi. Die Filterrichtung 

2 0 in a = - a max ist durch das Bild der im Pro j ektionswinkel 

a =— a max + Ls/2 zentrierten und in den Spiralgang eingepass- 
ten Segmentebene Ri gegeben. Fiir a > -a max + L s /2 ist das Bild 
A von der gekriimmten Flache Hi im virtuelle Detektor keine ge- 
jtifl rade Linie. In der Verof f ent lichung H. Turbell, P. E . 
^ ' Danielsson, An improved Pi-method for reconstruction from 

helical cone-beam projections, Proc . Int. Conf . On Fully 3D 
image Reconstruction, Egmond and Zee, The Netherlands, Jun 
23-26, 1999 wurde gezeigt, dass die Segmentebene Ri eine gute 
Naherung fur die gekriimmte Flache Hi darstellt . Daher wird 

30 sich in alien anderen Projektionen im Spiralsegment Ir 1 das 
Bild Hi im virtuellen Detektor urn die entsprechende Schnitt- 
linie von R x gruppieren. Dies bedeutet, dass fiir alle iibrigen 
Proj ektionswinkel im Teilsegment Ir 1 die Filterrichtungen von 
Di, die durch das Bild von Ri im virtuellen Detektor gegeben 

35 sind, geeignet sind. Obige Argumentation gilt gleichermafien 
fiir alle iibrigen Teilsegmente I R k (l<k< N t iit) . 
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Entsprechend diesem Grundgedanken schlagen die Erfinder vor, 
das an sich bekannte Verfahren zur Bilderstellung fur einen 
mehrzeiligen Spiral -Computertomographen zu verbessern. Be- 
kannt ist, dass : 

5 - ein Untersuchungsobj ekt durch ein rotierendes und sich in 
Richtung der Rotationsachse bewegendes Strahlenbiindel be- 
zuglich seines Absorptionsverhaltens spiralf ormig abgetas- 
tet und die gemessenen Absorptionsdaten durch einen fla- 
chig ausgebildeten Detektor gesammelt werden, 
10 - die gemessenen Absorptionsdaten auf einen virtuellen De- 
tektor, vorzugsweise im Drehzentrum der Rotation, proji- 
ziert und gefiltert werden, und 

zur Rekonstruktion des Absorptionswertes mindestens eines 
Voxels werden die gemessenen und gefilterten Daten verwen- 
15 det, die durch Strahlen erzeugt sind, welche dieses min- 

destens eine Voxel durchdringen . 

Die Verbesserung des Verfahrens liegt in der besonderen Aus- 
wahl der Filterrichtung, indem die Filterung der zur Rekon- 
20 struktion benutzten Daten im Bild des virtuellen Detektors in 
Richtung der Pro j ekt ion darauf abgebildeter Spiralsegmente, 
die durch die spiralf ormige Abtastung uber einen vorgegebenen 
Winkelbereich entstehen, erf olgt . 



Hierdurch ergibt sich, dass bei den Daten unterschiedlicher ^ J 

Rekonstruktionssegmente eines Voxels unterschiedliche und je- 
weils optimale Filterrichtungen auftreten. Wobei sich zusatz- 
lich der Vorteil ergibt, dass die Steigung der Spiale bezie- 
hungsweise die Wahl des Tischvorschubs beliebig ist. Es kann 
30 also in z-Richtung eine uberlappende Abtastung berucksichtigt 
werden. 



35 



In einer bevorzugte Ausfuhrung des Verfahrens findet die Fil- 
terung entlang der Schnittlinie von doppelt geneigten Ebenen 
im virtuellen Detektor statt . 
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Gemaft einer weiteren bevorzugten Ausfiihrung des Verfahrens 
betragt der vorgegebene Winkelbereich fur ein Spiralsegment 
der Lange L s ~ 240°. Wobei der Winkel 240° der Summe aus 
180°+2xFacherwinkel p max entspricht, der zur Rekonstruktion 
5 eines Bildvoxels erforderlich ist. 

Alternativ wird auch vorgeschlagen, vor der Filterung eine 
Parallelsortierung der Strahlen zur Bildung des virtuellen 
Detektors durchzufuhren, wobei hier der vorgegebene Winkelbe- 
10 reich fur ein Spiralsegment der Lange L s = 180° ist. 

Desweiteren konnen die durch die Spiralsegmente zumindest an- 
nahernd gebildeten Segmentebenen eine maximale Neigung auf- 
weisen, die durch die Bedingung begrenzt wird, dass an den 
beiden Enden des Spiralsegments innerhalb des Messfeldes 
Strahlen fur die Ebene vorhanden sein mussen. 

Zusatzlich ist es vorteilhaft, zur 3D-Riickpro j ektion ein Spi- 
ralsegment I x der Lange L I= [-a^x/ +ot max! mit ctmax =M * n/p in 
Ntiit iiberlappende Teilsegmente Ii k (l<k< N t ii t ) der Lange L s 
aquidistant einzuteilen, deren Schwerpunkte sich maximal um 
L s unterscheiden, wobei p dem eingestellten Pitch entspricht, 
so dass fur die entstehenden Subsegmente I R k (l<k< N t ii t ) 
gilt: 

I R k =I f k ;l<k<N tilt 

I R l = // U {- a V max -a max } 

~> rs T Ntilt T Ntilt , , f V \ 

3 0 I R = I j KJp max, a max) 

und im Rekonstrukt ionssegment I R k (l<k< N t iit) projektions- 
weise das zu einem Bildvoxel gehorige Pro j ektionsdatum im De- 
tektorbild D k ermittelt wird, wobei a^max den maximalen Winkel 
3 5 darstellt, den ein Strahl durch das Voxel V erreicht. Das 

Pro j ektionsdatum entspricht dem im Detektor gefundenen Pro- 
j ektionswert , der dem Abbild des Bildvoxels im Detektor ent- 
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spricht, wobei das Abbild der Schnittpunkt des Strahls vom 
Fokus durch das Bildvoxel im Detektor ist . 

Beriicksichtigt m'an zusatzlich den etwas unterschiedlich lan- 
gen Weg der Strahlen durch das Untersuchungsob j ekt in Abhan- 
gigkeit vom Konuswinkel des benutzten Strahlenbiindels , so er- 
weist es sich als vorteilhaft die gemessenen Absorptionsdaten 
in Abhangigkeit des Konuswinkels des erzeugenden Strahls in 
Richtung der Rotationsachse von Detektor und Strahlenquelle , 
vorzugsweise mit dem Kosinus ihres Konuswinkels, zu wichten. 

Der zur spiralf ormigen Abtastung des Untersuchungsob j ektes 
verwendete flachig ausgebildete Detektor kann mit einer Viel- 
zahl von Detektorelementen versehen sein, die matrixartig in 
Zeilen und Reihen angeordnet sind. Allerdings besteht grund- 
satzlich auch die Moglichkeit anders gestalteter Verteilungen 
der Detektorelemente liber den Detektor, beispielsweise in wa- 
benformiger Anordnung oder in Anordnung mit unterschiedlich 
breiten Detektorzeilen, die gegebenenf alls auch Detektoren 
unterschiedlicher Ausdehnung in Rotationsrichtung und Stei- 
gungsrichtung haben konnen. 

Neben dem erfindungsgemaBen Verfahren schlagen die Erfinder 
auch ein CT-Gerat zur Abtastung eines Untersuchungsob j ekt s 
vor, welches mit einem von mindestens einem Fokus ausgehenden 
Strahlenbiindel und mit einem flachig ausgebildeten Detektor- 
array mit einer Vielzahl von verteilten Detektorelementen zum 
Detektieren der Strahlen des Strahlenbiindels ausgestattet 
ist, wobei sich der mindestens eine Fokus relativ zu dem Un- 
tersuchungsob j ekt auf mindestens einer das Untersuchungsob- 
jekt umlaufenden Fokusbahn mit gegenuberliegendem Detektor- 
array bewegt und zumindest Mittel zur Sammlung von Detektor- 
daten, Filterung und 3D-Ruckpro j ektion vorgesehen sind, wobei 
die Mittel zur Filterung derart gestaltet sind, dass das Ver- 
fahren gemali einem der voranstehenden Verf ahrensanspruche 
durchgefiihrt wird. Die genannten Mittel zur Filterung konnen 
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zumindest teilweise durch Programme oder Programm-Module ver- 
wirklicht werden. 

Insgesamt wird durch die Verwendung des erf indungsgemafien 
5 Verfahrens und des erf indungsgemafien CT-Gerates erreicht, 

dass durch die besondere Auswahl der Filterrichtung eine ver- 
besserte Bildqualitat erzielt wird. 

Im folgenden wird die Erfindung anhand eines/der bevorzugten 
10 Ausf uhrungsbeispiele (s) mit Hilfe der Figuren/der Figur naher 
beschrieben. Es zeigen im einzelnen: 

Mehrzeilen-CT mit perspektivischer Darstellung 
der Abtasteinheit ; 

Langsschnitt durch das Mehrzeilen-CT aus Figur 1; 
Projektion von Segmentebenen in der paralleli- 
sierten Darstellung eines physikalischen Detek- 
tors Dp auf einem virtuellen Detektor D c mit Dar- 
stellung von drei Schnittvolumina und zugehorigen 
Spiral segment en ; 
4d: Darstellung der Projektion der Segmentebenen 
auf virtuellen Detektoren mit unterschiedlichem 
Pro j ektionswinkel ; 

Darstellung des Verlaufs von Projektionsebenen an 
Hand eines 16-zeiligen Detektors und Anderung der 
Filterrichtung; 

Darstellung der Wanderung der Projektion eines 
Voxels iiber den Detektor bei f ort schreitendem 
Pro j ektionswinkel ; 

Bildliche Darstellung der aneinandergereihten Re- 
konstruktionssegmente zur Reprojektion eines Vo- 
xels . 

Die Figuren 1 und 2 zeigen eine teils perspektivische Dar- 
35 stellung eines zur Durchfuhrung des erf indungsgemafien Verfah- 
rens geeigneten Mehr zeilen-CT-Gerates der dritten Generation. 
Die mit 1 bezeichnete Messanordnung (=Gantry) weist eine 




Figur 1 : 



15 Figur 2 : 
Figur 3 : 



20 

Figuren 4a- 
Figur 5 : 
Figur 6 : 
3 0 Figur 7: 
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Rontgenstrahlenquelle 2 mit einer dieser vorgelagerten quel- 
lennahen Strahlenblende 3 auf. Das flachig ausgebildete Array 
mit mehreren Zeilen und Spalten von Detektorelementen ausge- 
bildetes Detektorsystem 5 mit einer diesem vorgelagerten de- 
tektornahen Strahlenblende 6 ist in der Figur 2 im Schnitt 
dargestellt. In der Figur 1 sind zur besseren Ubersichtlich- 
keit nur acht Zeilen L x bis L Q von Detektorelementen 4 darge- 
stellt. Das Detektorsystem 5 kann jedoch auch eine andere, 
vorzugsweise groliere, Zeilenzahl aufweisen ohne den Rahmen 
der Erfindung zu verlassen. Ebenso ist auch eine andere fla- 
chige Anordnung der Detektoren moglich. 

Die Rontgenstrahlenquelle 2 mit der Strahlenblende 3 einer- 
seits und das Detektorsystem 5 mit der Strahlenblende 6 ande- 
rerseits sind an einem Drehrahmen 7 einander derart gegen- 
uberliegend angebracht, dass ein im Betrieb des CT-Gerats von 
der Rontgenstrahlenquelle 2 ausgehendes, durch die einstell- 
bare Strahlenblende 3 eingeblendetes , pyramidenf ormiges Ront- 
genstrahlenbiindel, dessen Randstrahlen mit 8 bezeichnet sind, 
auf das Detektorsystem 5 auftrifft. Dabei ist die Strahlen- 
blende 6 dem mittels der Strahlenblende 3 eingestellten Quer- 
schnitt des Rontgenstrahlenbiindels entsprechend so einge- 
stellt, dass entsprechend unterschiedlicher Betriebsmodi nur 
derjenige Bereich des Detektorsystems 5 freigegeben ist, der 
von dem Rontgenstrahlenbiindel unmittelbar getroffen wird. In 
den Figuren 1 und 2 werden nur acht Zeilen von Detektorele- 
menten 4, verwendet, die weiteren punktiert angedeuteten Zei- 
len sind von der Strahlenblende 6 abgedeckt und daher nicht 
aktiv. 

Das Rontgenstrahlenbiindel weist einen Conewinkel cp auf, bei 
dem es sich urn den Of f nungswinkel des Rontgenstrahlenbiindels 
in einer die Systemachse Z und den Fokus F enthaltenden Ebene 
handelt. Der Of f nungswinkel des Rontgenstrahlenbiindels in ei- 
ner rechtwinklig zu der Systemachse Z liegenden und den Fokus 
F enthaltenden Ebene ( Facheroff nungswinkel ) betragt 2/3tn ax . 
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Der Drehrahmen 7 kann mittels einer Antriebseinrichtung 22 urn 
eine mit Z bezeichnete Systemachse in Rotation versetzt wer- 
den. Die Systemachse Z verlauft parallel zu der z-Achse eines 
in der Figur 1 dargestellten und raumlich rechtwinkligen Ko- 
ordinatensystems . 

Die Spalten des Detektorsystems 5 verlaufen ebenfalls in 
Richtung der z-Achse, wahrend die Zeilen quer zu der System- 
achse Z verlaufen. 

Urn ein Untersuchungsobj ekt , z.B. einen Patienten, in den 
Strahlengang des Rontgenstrahlenbiindel bringen zu konnen, ist 
eine Lagerungsvorrichtung 9 vorgesehen, die parallel zu der 
Systemachse Z, also in Richtung der z-Achse, verschiebbar 
ist, und zwar derart, dass eine Synchronisation zwischen der 
Rotationsbewegung des Drehrahmens 7 und der Translationsbewe- 
gung der Lagerungsvorrichtung in dem Sinne vorliegt. Das Ver- 
haltnis von Translations- zu Rotationsgeschwindigkeit ist 
einstellbar, indem ein gewunschter Wert fur den Vorschub p 
(=pitch) der Lagerungsvorrichtung 9 pro Umdrehung des Dreh- 
rahmens gewahlt wird. 

Es kann also ein Volumen eines auf der Lagerungsvorrichtung 9 
befindlichen Untersuchungsobj ekts im Zuge einer Volumenabtas- 
tung untersucht werden, wobei die Volumenabtastung in Form 
einer Spiralabtastung in dem Sinne vorgenommen werden kann, 
dass unter gleichzeit iger Rotation der Messeinheit 1 und 
Translation der Lagerungsvorrichtung 9 mittels der Messein- 
heit pro Umlauf der Messeinheit 1 eine Vielzahl von Projekti- 
onen aus verschiedenen Pro j ekt ionsrichtungen auf genommen 
wird. Bei der Spiralabtastung bewegt sich der Fokus F der 
Rontgenstrahlenquelle relativ zu der Lagerungsvorrichtung 9 
auf einer Spiralbahn S um das Untersuchungsobj ekt . Um eine 
vollstandige Rekonstruktion eines CT-Bildes pro Zeile von De- 
tektorelementen zu gestatten, muss die Spiralabtastung sich 
in a-Richtung liber wenigstens 7r+2/ibax erstrecken, sie kann in- 
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nerhalb cler technischen Grenzen des CT-Gerats aber auch be- 
liebig langer sein. 

Ein Volumen des Untersuchungsobj ekts kann jedoch infolge des 
Umstandes, dass mehrere Zeilen von Detektorelementen 4 vor- 
handen sind, auch im Zuge einer sogenannten Tomogrammab- 
tastung untersucht werden, bei der keine Relativbewegung in 
Richtung der z-Achse zwischen Messeinheit 1 und Lagerungsvor- 
richtung 9 stattfindet (p=0) . Im Falle der Tomogrammabtastung 
ist also die Grofte des untersuchten Volumens durch die Zahl 
der aktiven Zeilen von Detektorelementen 4 bestimmt. Wahrend 
einer Tomogrammabtastung bewegt sich der Fokus F auf einer 
kreisf ormigen oder bei vorhandenem Vorschub spiralf ormigen 
Fokusbahn, die in einer im Folgenden als Mittelebene bezeich- 
neten Ebene liegt. 

Die Tomogrammabtastung kann in Form eines Teilumlaufs oder in 
Form eines Vollumlaufs erfolgen, wobei der Teilumlauf ein 
Teilumlauf intervall von wenigstens ^"+2>ft l ax (ein Halbumlauf 
plus Facherof f nungswinkel) umfasst, wodurch eine vollstandige 
Rekonstruktion eines CT-Bildes gestattet wird, wahrend ein 
Vollumlauf 2k umf asst . 

Die wahrend der Spiral- oder Tomogrammabtastung aus den De- 
tektorelementen jeder aktiven Zeile des Detektorsystems 5 pa- 
rallel ausgelesenen, den einzelnen Projektionen P(a,/3,q) in 

Facherstrahlgeometrie entsprechenden Messdaten werden in ei- 
ner Datenauf bereitungseinheit 10 einer Digital/Analog-Wand- 
lung unterzogen, serialisiert und an einen Bildrechner 11 
ubertragen . 

Nach einer Vorverarbeitung der Messdaten in einer Vorverar- 
beitungseinheit 12 des Bildrechners 11 gelangt der resultie- 
rende Datenstrom zu einer Schnittbildrekonstrukt ionseinheit 
13, die aus den Messdaten Schnittbilder von gewunschten 
Schichten des Untersuchungsobj ekts in an sich bekannter Weise 
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nach einem 3D~Ruckpro j ektionsverf ahren rekonstruiert, wobei 
jedoch die verwendeten Daten zunachst gefiltert werden. 

Die CT-Bilder setzen sich aus matrixartig zusammengeset zten 
5 Pixeln (Pixel= picture element) zusammen, wobei die Pixel der 
jeweiligen Bildebene zugeordnet sind, jedem Pixel eine CT- 
Zahl in Hounsfield Units (HU) zugeordnet ist und die einzel- 
nen Pixel entsprechend einer CT-Zahl/Grauwertskala in einem 
ihrer jeweiligen CT-Zahl entsprechenden Grauwert dargestellt 
10 werden. Dabei veranschaulicht jedes Pixel ein Voxel (Voxel = 
volume element ) der in dem CT-Bild veranschaulichten Schicht 
des Untersuchungsobj ekts . Da infolge der Mehrzeiligkeit des 
■ Detektorsystems 5 und gegebenenf alls der Spiralabtastung 

Messdaten beziiglich mehrerer Schichten des Untersuchungsob- 
15 jekts gewonnen werden, stehen 3D-Daten zur Verfligung, die im 
Rahmen der Erfindung einer 3D-Riickpro j ektion unterzogen wer- 
den. Als Endergebnis stehen 3D-Bilddaten in Form einer drei- 
dimensionalen Matrix, beispielsweise mit den Achsen x, y, z 
zur Verfugung, wobei jedes Element der Matrix einem Voxel 
20 V{ X ,y,z) entspricht und den der zugehorigen CT-Zahl entspre- 
chenden Grauwert enthalt. Diejenigen Elemente der dreidimen- 
sionalen Matrix, die den gleichen x- , y- , oder z-Wert aufwei- 
sen, stellen dann jeweils ein planares Schnittbild der dem 
maftgeblichen x- , y- , oder z-Wert entsprechenden Schicht des 
Untersuchungsobj ekts dar . 

Die von der Schnittbildrekonstruktionseinheit 13 rekon- 
struierten Bilder werden auf einer an den Bildrechner 11 an- 
geschlossenen Anzeigeeinheit 16, z.B. einem Monitor, darge- 
30 stellt. 

Die Rontgenstrahlenquelle 2, beispielsweise eine Rontgen- 
rohre, wird von einer Generatoreinheit 17 mit den notwendigen 
Spannungen und Stromen, beispielsweise der Rohrenspannung U, 
35 versorgt. Urn diese auf die jeweils notwendigen Werte einstel- 
len zu konnen, ist der Generatoreinheit 17 eine Steuereinheit 
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18 mit Tastatur 19 zugeordnet, die die notwendigen Einstel- 
lungen gestattet . 

Auch die sonstige Bedienung und Steuerung des CT-Gerates er- 
folgt mittels der Steuereinheit 18 und der Tastatur 19, was 
dadurch veranschaulicht ist, dass die Steuereinheit 18 mit 
dem Bildrechner 11 verbunden ist . 

Unter anderem kann die Anzahl der aktiven Zeilen von Detek- 
torelementen 4 und damit die Position der Strahlenblenden 3 
und 6 eingestellt werden, wozu die Steuereinheit 18 mit den 
Strahlenblenden 3 und 6 zugeordneten Verstelleinheiten 2 0 und 
21 verbunden ist. Weiter kann die Rotat ionszeit r eingestellt 
werden, die der Drehrahmen 7 fur eine vollstandige Umdrehung 
benotigt, was dadurch veranschaulicht ist, dass die dem Dreh- 
rahmen 7 zugeordnete Antriebseinheit 22 mit der Steuereinheit 
18 verbunden ist. 

Obwohl es grundsat zlich auch moglich ist, das erf indungsge- 
mafie Verfahren auch in Facherstrahlgeometrie zu realisieren, 
wird das beschriebene CT-Gerat vorzugsweise in einem Modus 
betrieben, in dem das erf indungsgemafie Verfahren in Parallel- 
strahlgeometrie realisiert ist. 

Demnach werden die bei der Abtastung des fur die jeweilige 
Untersuchung relevanten Korperbereichs des Patienten durch 
Spiral- oder Tomogrammabtastung in Facherstrahlgeometrie ge- 
wonnenen Daten zunachst in an sich bekannter Weise durch ein 
im Allgemeinen als " rebinning" bezeichnetes Verfahren in Da- 
ten in Parallelstrahlgeometrie umgewandelt. Diese Umwandlung 
beruht auf einer Umsortierung der in Facherstrahlgeometrie 
gewonnenen Daten derart, dass aus unterschiedlichen in Fa- 
cherstrahlgeometrie auf genommenen Projektionen Strahlen ent- 
nommen und zu einer Projektion in Parallelstrahlgeometrie zu- 
sammengefugt werden. In Parallelstrahlgeometrie reichen Daten 
aus einem Intervall der Lange n aus, urn ein vollstandiges 
Bild rekonstruieren zu konnen. Urn diese Daten gewinnen zu 
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konnen, mussen Daten in Facherstrahlgeometrie aus einem In- 
tervall der Lange ;r+2y#n a x zur Verfiigung stehen. 

Die Figur 3 zeigt die Spiralbahn S der Messanordnung des Spi- 
ral-CT, auf der sich der Fokus F mit dem gegenuberliegenden 
physikalischen Detektor D p bewegt , wobei in dieser Darstel- 
lung bereits die Umwandlung in die Parallelstrahlgeometrie 
durchgefiihrt wurde. Erkennbar ist dies auch dadurch, dass in 
diesem Bild der Fokus F in ftinf nebeneinander liegenden Posi- 
tionen mit je drei parallelen Strahlen gezeigt ist. Diese Pa- 
rallelstrahlen treffen im Hintergrund auf den physikalischen 
Detektor D p/ der allerdings auch schon lediglich ein Abbild 
von vielen Detektorpositionen umlaufender Detektorelemente 
darstellt. Das Detektorbild des physikalischen Detektors D p 
ist weiterhin auf den virtuellen Detektor D v projiziert, der 
an der Systemachse anliegt. Aulierdem zeigt die Figur 3 drei 
beispielhaf ten Schnittebenen E x bis E 3 , die auf dem Detektor 
D v abgebildet werden. 

Betrachtet man nun den virtuellen Detektor D v im Koordinaten- 
system alleine mit der Darstellung der Projektionen P x bis 
P 3/ die den Ebenen Ei bis E 3 aus Figur 3 entsprechen, so er- 
halt man eine Darstellung des virtuellen Detektors D v , wie 
sie in den Figuren 4a bis 4d gezeigt ist. Die sukzesive Ver- 
schiebung der Projektionen in den Figuren 4a bis 4d entsteht 
dabei durch die Betrachtung der Ebenen El bis E3 und deren 
Projektionen PI bis P3 bei f ortlauf ender Drehung um die z- 
Achse . 

Der erf indungsgemafte Gedanke besteht nun darin, auch in der 
3D-Rekonstruktion eine Filterung durchzuf iihren, die in der 
jeweiligen Pro j ektionsrichtung der betrachteten Ebene ver- 
lauft. Wie aus den Figuren 4a bis 4d ersichtlich ist, andert 
sich die Ausrichtung der Projektionen von unten nach oben 
kontinuierlich und somit ebenfalls die Filterrichtung . 
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Betrachtet man eine entsprechende Darstellung mit einer hohe 
ren Anzahl von Ebenen und Pro j ektionen, so erhalt man eine 
Darstellung des virtuellen Detektors, wie sie in Figur 5 ge- 
zeigt ist. Die Pro j ektionen P andern dabei ihre Richtung von 
unten nach oben in typischer kontinuierlicher Weise. Entspre 
chend wird erf indungsgemaft fur jedes so gewonnene Detektor- 
bild eine zweidimensionale Filterung von links nach rechts 
durchgef uhrt , deren Richtung sich auch uber die Hohe des De- 
tektor entsprechend Verlauf srichtung der Projektionen veran- 
dert . 

Zur 3D-Ruckpro j ektion, die sich auf die Voxel V (X/y , 2) bezieht, 
werden zur Bestimmung des Absorptionsvermogens eines Voxels 
die Strahlen betrachtet, die dieses Voxel im Untersuchungsvo 
lumen durchdr ingen . Beim kontinuierlichen Umlauf um die z- 
Achse und Projektionen dieses einen Voxels V (x , y/Z) auf den 
virtuellen Detektor D v verlauft die Bahn der Projektion die- 
ses Voxels uber die Vielzahl der betrachteten Detektorbilder 
wie es in der Figur 6 mit der Bahn B v dargestellt ist. Hier 
sind fur unterschiedliche Pro j ektionsrichtungen die Schnitt- 
linie einer Segmentebene im jeweiligen virtuellen Detektor 
eingetragen. In diesem Fall werden 32 unterschiedliche Pro- 
jektionsrichtungen betrachtet. Durch das erf indungsgemafte 
Verfahren erhalten die einzelnen Detektorpunkte dabei eine 
Filterung in unterschiedliche Richtungen f, die jeweils der 
Richtung der Projektionen der Segmentebene, welche durch die 
ses Voxel verlauft, entspricht. Zur Vermeidung von Missver- 
standnissen ist anzumerken, dass diese Anzahl nicht mit der 
Zeilenzahl des Detektors gleichzuset zen ist. 

Die dargestellte Bahn B v des Voxels entspricht einem Umlauf 
der Messanordnung zwischen -amax V und +a max v - Teilt man diese 
Bahn in erf indungsgemafter Weise in Rekonstruktionsunterseg- 
mente Ir 1 bis I R Ntllt auf, so lassen sich diese wie in Figur 7 
gezeigt darstellen. Hierbei ist darauf hinzuweisen, dass der 
dargestellte Gesamtwinkel von -a raax v und +a max v natiirlich weit 
grolier als die gezeigten 180° ist und mehreren Umlauf en der 
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Messanordnung entspricht. Die Segmente sind dabei - mit Aus- 
nahme der randstandigen Segmente - bezuglich ihrer Mittel- 
und Schwerpunkte aquidistant und leicht uberlappend angeord- 
net . 



Die Riickpro j ektion dieser Rekonstruktionssegmente auf das be- 
trachtete Voxel V (X/y/Z) ergibt dann den diesem Voxel zugeord- 
neten Absorptionswert innerhalb des Untersuchungsobj ektes . 
Wird diese Vorgehensweise fur alle Voxel im gesamten betrach- 
10 teten Untersuchungsvolumen durchgef iihrt , so erhalt man eine 
Volumendarstellung der Absorptionswerte des Untersuchungsob- 

* 3ektes - 

Weil jedes Bildvoxel in jedem Teilsegment des Rekonstrukt i - 

• 15 onsbereiches im Detektor abgebildet wird, tragt jedes Bildvo- 

xel mit dem gleichen Gewicht zum rekonstruierten Bildvolumen 
bei. Daher ist eine voxelweise Gewichtung der gemessenen Da- 
ten nicht notwendig. 

20 Insgesamt wird also mit dieser Erfindung ein Verfahren zur 



Bilderstellung fur einen mehrzeiligen Spiral-CT mit 3D-Riick- 
proj ektion und ein Spiral-CT beschrieben, bei dem zur Rekon- 
struktion des Absorptionswertes mindestens eines Voxels die 
gemessenen und gefilterten Daten verwendet werden, die durch 



5 Strahlen erzeugt sind, welche dieses mindestens eine Voxel 



durchdringen, wobei die Filterung der zur Rekonstruktion be- 
nutzten Daten im Bild des virtuellen Detektors in Richtung 
der Proj ektion darauf abgebildeter Spiralsegmente , die durch 
die spiralf ormige Abtastung uber einen vorgegebenen Winkelbe- 
30 reich entstehen, erf olgt . 

Es ist erganzend darauf hinzuweisen, dass das zuvor beschrie- 
bene Verfahren und CT auch zur Herstellung von periodisch be- 
wegten Objekten, insbesondere zur Cardio-Spiral-CT, verwendet 
35 werden kann, wobei in diesem Falle in an sich bekannter Weise 
eine Gewichtung oder Auswahl relevanter Daten in Korrelation 



5 
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zur Bewegungssituation des Untersuchungsob j ektes durchgef iihrt 
werden muss. 
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Patentanspriiche 

1. Verfahren zur Bilderstellung fur einen mehr zeiligen 
Spiral - Compu t e r t omogr aphen , wobei 

1.1. ein Untersuchungsobj ekt durch ein rotierendes und sich 
in Richtung der Rotat ionsachse bewegendes Strahlenbiin- 
del bezuglich seines Absorptionsverhaltens spiralf ormig 
abgetastet und die gemessenen Absorptionsdaten durch 
einen flachig ausgebildeten Detektor gesammelt werden, 

1.2. die gemessenen Absorptionsdaten auf einen virtuellen 
Detektor, vorzugsweise im Drehzentrum der Rotation, 
projiziert und gefiltert werden, und 

1.3. zur Rekonstrukt ion des Absorptionswertes mindestens ei- 
nes Voxels werden die gemessenen und gefilterten Daten 
verwendet, die durch Strahlen erzeugt sind, welche die- 
ses mindestens eine Voxel durchdringen, 

dadurch gekennzeichnet, dass 

1.4. die Filterung der zur Rekonstruktion benutzten Daten im 
Bild des virtuellen Detektors in Richtung der Projek- 
tion darauf abgebildeter Spiralsegmente , die durch die 
spiralf ormige Abtastung iiber einen vorgegebenen Winkel- 
bereich entstehen, erf olgt . 

2. Verfahren gemaft dem voranstehenden Patentanspruch 1, 
dadurch g e ke nn z e i c hne t , dass die Filterung ent- 
lang der Schnittlinie von doppelt geneigten Ebenen im 
virtuellen Detektor stattfindet. 

3 . Verfahren gemaft einem der voranstehenden Patentansprii- 
che 1 bis 2, dadurch gekennzeichnet , dass der 
vorgegebene Winkelbereich fiir ein Spiralsegment der 
Lange LS < n+2*p max betragt. 

4 . Verfahren gemafi einem der voranstehenden Patentansprii- 
che 1 bis 2, dadurch gekennzeichnet, dass vor 
der Filterung eine Parallelsortierung der Strahlen zur 
Bildung des virtuellen Detektors stattfindet. 
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5. Verfahren gemafi dem voranstehenden Patentanspruch 4, 
dadurch gekennzeichnet , dass der vorgegebene 
Winkelbereich fur ein Spiralsegment der Lange LS < 180° 

5 betragt. 

6. Verfahren gemafi einem der voranstehenden Patentansprii- 
che 1 bis 5, dadurch gekennzeichnet, dass die 
durch die Spiralsegmente zumindest annahernd gebildeten 

10 Segmentebenen eine maximale Neigung besitzen, so dass 

an den Enden des betrachteten Spiralsegments innerhalb 
des Messfeldes Strahlen fur die Segmentebene im Detek- 
tor vorhanden sind. 

-15 7. Verfahren gemaft einem der voranstehenden Patentansprii- 
che 1 bis 5, dadurch gekennzeichnet, dass zur 
3D-Ruckprojektion ein Spiralsegment I T der Lange 
Li= [-octnax/ +oc max ] mit Otnax = M * n /p in N t ii t uberlappende 
Teilsegmente (l<k< N ti i t ) der Lange L s aquidistant 

20 eingeteilt wird, deren Schwerpunkte sich maximal um L s 

unterscheiden, wobei p dem eingestellten Pitch ent- 
spricht, so dass fur die entstehenden Subsegmente I R k 
(l<k< Ntiit) gilt: 

^25 I R k =I f k ;l<k<N tilt 



30 



I R — I j W |— Ct V imx 9 — CC max} 

T Ntilt j Ntilt ( V \ 

I R = Ij )Ct max, a max) 



und im Rekonstruktionssegment I R k (l<k< N t iit) projekti- 
onsweise das zu einem Bildvoxel gehorige Pro j ektionsda- 
tum im Detektorbild D K ermittelt wird, wobei a^max den 
maximalen Winkel darstellt, den ein Strahl durch das 
35 Voxel V erreicht. 
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8. Verfahren gemaft einem der voranstehenden Patentanspru- 
che 1 bis 7, dadurch gekennzeichnet , dass die 
gemessenen Absorptionsdaten in Abhangigkeit des Konus- 
winkels des erzeugenden Strahls in Richtung der Rotati- 

5 onsachse von Detektor und Strahlenquelle , vorzugsweise 

mit dem Kosinus ihres Konuswinkels , gewichtet werden. 

9. Verfahren gemaft einem der voranstehenden Patentansprii- 
che 1 bis 8, dadurch gekennzeichnet, dass zur 

10 spiralf ormigen Abtastung ein flachig ausgebildeter De- 

tektor mit einer Vielzahl von Detektorelementen verwen- 
det wird, die matrixartig in Zeilen und Reihen angeord- 
net sind. 

15 10. CT-Gerat zur Abtastung eines Untersuchungsob j ekt s mit 
einem von mindestens einem Fokus ausgehenden Strahlen- 
bundel und mit einem flachig ausgebildeten Detektorar- 
ray mit einer Vielzahl von verteilten Detektorelementen 
zum Detektieren der Strahlen des Strahlenbiindels , wobei 
20 der mindestens eine Fokus relativ zu dem Untersuchungs- 

objekt auf mindestens einer das Untersuchungsob j ekt urn- 
laufenden Fokusbahn mit gegeniiberliegendem Detektor- 
array bewegt, wobei zumindest Mittel zur Sammlung von 
Detektordaten, Filterung und Ruckpro j ektion vorgesehen 
sind und die Mittel zur Filterung derart gestaltet 
sind, dass das Verfahren gemafi einem der voranstehenden 
Verf ahrensanspruche durchgef iihrt wird . 

11. CT-Gerat gemaft dem vorangehenden Anspruch, dadurch 
30 gekennzeichnet, dass die genannten Mittel zur 

Filterung zumindest teilweise durch Programme oder Pro- 
gramm-Module verwirklicht werden. 



200210626 ' , 

21 

Zusammenf as sung 

Verfahren zur Bilderstellung fur einen mehrzeiligen Spiral- 
Computertomographen mit 3d-Rekonstruktion und Computer- 
Tomographie-Gerat zur Durchfuhrung dieses Verfahrens 

Die Erfindung betrifft ein Verfahren zur Bilderstellung fur 
einen mehrzeiligen Spiral-CT mit 3D-Ruckproj ektion, wobei zur 
Rekonstruktion des Absorptionswertes mindestens eines Voxels 
die gemessenen und gefilterten Daten verwendet werden, die 
durch Strahlen erzeugt sind, welche dieses mindestens eine 
Voxel durchdringen . Aufierdem betrifft die Erfindung ein CT- 
Gerat zur Durchfuhrung dieses Verfahrens mit entsprechenden 
Mitteln zur Bilderstellung. 

Beim erf indungsgemafien Verfahren erfolgt die Filterung der 
zur Rekonstruktion benutzten Daten im Bild des virtuellen De- 
tektors in Richtung der Pro j ektion darauf abgebildeter Spi- 
ralsegmente, die durch die spiralf ormige Abtastung uber einen 
vorgegebenen Winkelbereich entstehen. Das erf indungsgemafie 
CT-Gerat weist Mittel auf, eine gerichtete Filterung dieser 
Art durchfiihren. 



Figur 3 
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